Global asymmetry of many-qubit correlations: A lattice gauge theory approach 
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We introduce a novel bridge between the familiar gauge field theory approaches used in many ar- 
eas of modern physics such as quantum field theory and the SLOCC protocols familiar in quantum 
information. Although the mathematical methods are the same the meaning of the gauge group will 
be difi^erent. The measure we introduce, 'twist', is constructed as a Wilson loop from a correlation 
induced holonomy. The measure can be understood as the global asymmetry of the bipartite corre- 
lations in a loop of three or more qubits; if the holonomy is trivial (the identity matrix), the bipartite 
correlations can be globally untwisted using general local qubit operations, the gauge group of our 
theory, which turns out to be the group of Lorentz transformations familiar from special relativ- 
ity. If it is not possible to globally untwist the bipartite correlations in a state globally using local 
operations, the twistedness is given by a non-trivial element of the Lorentz group, the correlation 
induced holonomy. We provide several analytical examples of twisted and untwisted states for three 
qubits, the most elementary non-trivial loop one can imagine. 

PACS numbers: 03.67.Mn, 03.65.Vf, ll.15.Ha, ll.30.Cp 



I. INTRODUCTION 



There has been a large effort during the last decade 
to study the properties of entanglement and correlations 
in quantum states. Progress has been made quantify- 
ing and characterizing these correlations for composite 
states of two and more qubits [l|. However, for states 
of three or more qubits the problem becomes hard as 
the possible ways a state may be entangled becomes 
greater. In this paper we introduce a new approach in 
the spirit of a lattice gauge field theory. Our motivation 
to use a gauge theory approach comes from several lines: 
(i) Our most successful theories of nature are expressed 
in this language and are well developed mathematically. 
Some examples of gauge theories in physics include clas- 
sical electromagnetism, quantum field theories (such as 
quantum electrodynamics and quantum chromodynam- 
ics), string theory and general relativity. Mapping the 
study of quantum correlations to a gauge theory may 
allow one to use some of the techniques previously de- 
veloped, (ii) Gauge theories have a natural geometric 
interpretation with which to gain intuition about a prob- 
lem, the gauge field defining a curved surface and the 
Wilson loop, a gauge invariant observable of the theory, 
giving a measure of the total curvature of the gauge field. 
In our work, the gauge group that emerges naturally is 
the Lorentz group. 

The usual scenario one considers when studying entan- 
glement is the following: There are N parties in spatially 
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separated locations each holding one part of a composite 
quantum state comprised of N subsystems. Each party 
is free to make local operations on their part of the state 
and classically communicate this operation and its out- 
come to the other parties in an effort to put the state in 
a standard form. They can use deterministic operations 
given by unitary groups which are simple rotations of 
their subsystems or more general local operations given 
by Kraus operators that may only be successful proba- 
bilistically. If two different states can be brought to the 
same standard form reversibly they are said to belong to 
the same entanglement class. If the conversion is deter- 
ministic, the two states are unitarily equivalent. If the 
conversion is only probabilistic, the two states are said 
to be equivalent under stochastic local operations and 
classical communication (SLOCC) 0- In this paper we 
concentrate on the latter scenario for qubits. In this case 
the group of reversible local qubit operations is SL{2, C) 
up to a positive constant less than or equal to unity. 
Moreover, it is known that the entanglement measures 
concurrence W\ and three-tangle 14] are invariant under 
the action of S'i (2, C). 

There are a number of results concerning SLOCC clas- 
sification. Diir et al. ^ have shown that there exist two 
classes for pure states of three qubits, the Greenberger- 
Horne-Zeilinger (GHZ) class and the W class. Verstraete 
et al. have shown that there exist nine families for 
pure states of four qubits (although there are infinite 
number of classes; one of these families depends on a 
continuous parameter). If one is given a single copy of 
a state and allowed to use SLOCC, Lo and Popescu 
have shown that any pure two qubit state can be brought 
to a maximally entangled Bell state and Kent et al. and 
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others [8l-ll2l| have shown that the most entangled state 
one can obtain from a mixed state of two qubits is Bell 
diagonal. There has also been work on creating entangle- 
ment monotones by exploiting the SL(2, C) invariance of 
certain entanglement measures |13l4l5| . 

In a lattice gauge theory the gauge field assigns a trans- 
formation to every pair of neighboring lattice points [l6l | . 
This transformation is an element of the gauge group 
and is sometimes known as a parallel transporter; it par- 
allel transports the property sitting at one lattice point 
to the neighboring point. For example, in lattice quan- 
tum chromo dynamics, space-time is discretized into lat- 
tice points and the property sitting on the lattice points 
is color represented by a vector. To see how color changes 
from one lattice point to another one multiplies the vec- 
tor by the parallel transporter. The assignment of paral- 
lel transporters to each link specifies the configuration of 
the gauge field. While the parallel transporters them- 
selves change under gauge transformations, there is a 
natural gauge invariant observable on the lattice given 
by the trace of the total transformation around a loop 
(or plaquette). This is known as the Wilson loop. This 
total transformation is a measure of the curvature of 
the gauge field around that loop (equivalently the flux 
through the loop). The closer to the identity the trans- 
formation around the loop is, the less curved the gauge 
field is. 

In this paper our N lattice points are the N qubits 
of a composite quantum state, the transformations be- 
tween neighboring qubits are specified by the local oper- 
ations that symmetrize the correlations in that two-qubit 
link and the gauge invariant observable, the Wilson loop, 
provides a measure of the degree of asymmetry of the bi- 
partite correlations globally. We call this measure twist. 
These are correlations such as entanglement invariant un- 
der SL{2, C), the gauge group of our theory. In an earlier 
paper, Wootters explored a related idea and found evi- 
dence that a non-trivial twisting around a loop requires 
the sacrifice of some entanglement. The definitions of 
twist that we introduce in the present paper are different 
from the one defined in the earlier paper and they 
do not have the same interpretation. 

The interpretation of our measure is as follows: Imag- 
ine taking a state comprised of N qubits. One performs a 
local Kraus operation on the second qubit to symmetrize 
the correlations between the first and second qubits (sym- 
metrize here means that the correlations specified by the 
expectation values of spin measurements on the first and 
second qubits labeled 1 and 2 are the same under inter- 
change of the two qubits i.e. {aj cr|) = {aj af)). The 
indices i and j take the values 0, 1, 2, 3 with ctq = I and 
<7i, <T2, CT3 are the Pauli matrices. We can now imagine 
performing a local operation on the third qubit to sym- 
metrize the correlations between qubits two and three, 
the link 2, 3. Now the links comprised of qubits 1, 2 and 
2, 3 have been symmetrized. We can repeat this pro- 
cess along all links until the final one comprised of the 
last qubit, qubit N and the first qubit. One now has 



a dilemma; the local operation performed on the first 
qubit that symmetrizes the link N, 1 may cause the link 
I, 2 to become asymmetric. The degree of mismatch be- 
tween the initial and final local operations is the total 
transformation (also known as a holonomy) of an under- 
lying correlation-induced gauge field around the loop and 
a measure of this mismatch is given by the trace of this 
total transformation, the Wilson loop. If one can simul- 
taneously symmetrize all two qubit links in the loop, the 
overall transformation is trivial i.e. the identity. How- 
ever, if this is not possible, the Wilson loop gives a degree 
of asymmetry in the bipartite correlations. 

Even though our measure is defined for states of N 
qubits, in this paper we restrict study to three qubit 
states, the most elementary loop one can conceive of. 
Explicitly, we show that all two qubit states have a triv- 
ial total transformation whereas all pure states of three 
qubits have a total transformation equal to a tt rotation. 
We then look at several examples of mixed states of three 
qubits. We provide examples of untwisted states and two 
examples of states with 5*0(1, 1) holonomy and link them 
with the concurrence in each link. These are examples we 
have been able to work out analytically since 50(1, 1) is 
a simple group depending on one parameter; it appears 
however, a generic mixed state of three qubits has the 
full 50+(l, 3) structure depending on all six parameters. 
That is, a generic state cannot be untwisted. Although 
some analytical examples may be laborious to work out, 
one can calculate twist easily and quickly using numerical 
methods. 

The ideas we outline here provide a novel bridge be- 
tween the familiar gauge field theory approaches used 
in many areas of modern physics such as quantum field 
theory and the SLOCC protocols familiar in quantum 
information. The mathematical methods are the same 
although the meaning is very different. For example, the 
geometric phase [l8| has the same mathematical struc- 
ture as our measure twist, however the holonomy is a 
result of the curvature of the space of quantum states 
or the parameter space of a Hamiltonian. In the case of 
twist the holonomy is induced by the bipartite correla- 
tions between quantum states. The key in the mapping 
of correlations to parallel transporters in this paper is the 
introduction of a Lorentz polar decomposition. 

The paper is organized as follows: In the following sec- 
tion we show how one may assign a parallel transporter 
to a two qubit link and we introduce the measure twist. 
We then provide analytical results for states of two and 
three qubit states in sections IIIII and IIVI respectively be- 
fore concluding in section [V] 



II. MAPPING THE CORRELATIONS IN A 
QUANTUM STATE TO A PARALLEL 
TRANSPORTER 

In this section we show how we map a parallel trans- 
porter to each two qubit link, the properties of these par- 
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allel transporters and the gauge invariance of a measure 
we caU 'twist', given by the Wilson loop. 

The approach we use here is inspired by our previous 
work |l9j which showed how to generate local invariants. 
We represented bipartite states by correlation matrices, 
sometimes known as the Hilbert-Schmidt representation. 
That is, the two qubit density matrix, pab is now repre- 
sented by a real (although generally not positive or sym- 
metric) four by four correlation matrix, S{a, b) whose 
elements i,j are given by 

5(a,% = itr[«®a^^)p,b]. (1) 

We can also imagine performing local operations on the 
state Pab to take it to a new state pj^^ = A(^BpabA^ (^B^ . 
In the Hilbert-Schmidt representation, the local opera- 
tions act on S{a, b) as follows 

S{a,by = AS{a,b)B^ (2) 

where the elements of the real, four by four matrices of 
the local operations on a and b in the new representation 
are given by 

A.,., = itr(^t^.^^a)^ 

eL = it^(5-).^MJ- (3) 

The local operation A must obey the constraint AA^ < I 
(and likewise for B) to take p to a state with normaliza- 
tion < 1. However, the standard practice when thinking 
about SLOCC protocols is to rescale A and i? by a posi- 
tive constant so they have determinant one. That is, they 
are now elements of SL{2, C) if they are reversible opera- 
tions. In a SLOCC protocol one does not care about nor- 
malization, only that it may be possible to convert one 
state to another with some non-zero probability. This 
rescaling also makes sense when studying entanglement 
as the concurrence and three-tangle in a state remain in- 
variant when considering local operations in this group as 
previously mentioned. From this point on we will restrict 
the local operations to be elements of SL{2,C). When 
making this rescaling, the local operations in the Hilbert- 
Schmidt basis now become elements of SO'^{l,3), the 
gro up of proper, orthochronous Lorentz transformations 
[111 |20| . This is due to the well known homomorphism 
SL{2,C) ~ S'0+(l,3). One can verify that A and B 
are indeed elements of V & S'0"'"(l,3) from the defin- 
ing property of the Lorentz group; that it preserves the 
Minkowski metric, rj = diag{l, —1, —1, —1}, 

Vr]V'^ = rj. (4) 

77 is equivalent to the spin-flip operator used in calcu- 
lation of the well known entanglement measure, concur- 
rence and is also equivalent to the singlet state |01) — |10) 
when written as a correlation matrix. 



A. Lorentz polar decomposition and association of 
a parallel transporter to each two-qubit link 

Given a single copy of an arbitrary two-qubit state 
what is the most entangled state one can convert it to 
using SLOCC? In a series of papers (8l-[T2j it was shown 
that the most entangled state one can produce is a Bell 
diagonal state, a mixture of some or all of the four Bell 
states = |01) ± |10), |$±) = |00) ± |11). In the 

correlation matrix representation, the four Bell states are 
the diagonal matrices 

-^1*-) = 5 diag{l, -1,-1,-1}, 

5*1*+) = \ diag{l, 1,1,-1}, 

5|<i>-) = \ diag{l, -1,1,1}, 

= i diag{l, 1,-1,1}. (5) 

Since a Bell diagonal state is just a mixture of these 
pure Bell state correlation matrices, a Bell diagonal state 
must also be a diagonal correlation matrix. This prompts 
a natural decomposition of an arbitrary two-oubit state 
represented by a correlation matrix S{a,b) [ll| : 

Sia,b)^Va^abW;^, (6) 

where = diagjso, si, S2, S3} represents a Bell diag- 
onal state and Va and Wb (elements of S'O''" (1,3)) are 
the local operations taking one from the Bell diagonal 
state to the state represented by S{a,b). This is a spe- 
cial form of a singular value decomposition (SVD) with 
a SLOCC operational interpretation. With the usual 
SVD one makes the entries of S, the singular values, 
real and positive and puts them in non-increasing order 
down the diagonal. Since V,W ^ S'0+ (1,3) we can or- 
der the Lorentz singular values in non-increasing order 
and they are real since S, V and W are real. However, 
we cannot generally make them all positive since doing 
so could result in a density matrix with negative prob- 
abilities, an invalid physical state. The leading Lorentz 
singular value, sq, is always positive and the largest (it 
represents the normalization of the state) . The remain- 
ing three Lorentz singular values si, S2 and S3 take the 
same sign (generically negative). The ordering of these 
remaining three singular values will not affect the paral- 
lel transporter we assign to each link, however the sign 
they take will. These singular values represent the ex- 
pectation values of spin in the three spatial directions 
Si = so{ai <Ji) rescaled by the normalization of the 
state Sq. 

The Lorentz singular values are SL(2, C) invariants 
and are closely related to the concurrence in the two- 
qubit state [10, El. That is, the entanglement in the 
state is related to the (dilated) magnitude of the corre- 
lations. It is a Lorentz scalar. The concurrence C for an 
arbitrary two qubit state is 

C = max{0, -trS}. (7) 
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By rewriting the Lorentz SVD, we can also find a left 
Lorentz polar decomposition. That is, we may re-write 
S as 



5(a,6) = A(a,&)5(a,6) (8) 



where A(a,6) = VavWjfr) e SO+{l,3) and S{a,b) = 
WbTiabWif . This is equivalent to our previous decom- 
position of S except now S{a, b) represents a state with 
correlations that are symmetric under interchange of the 
qubits and A(a, b) represents the local operation per- 
formed on a that symmetrizes the correlations. More- 
over if S is full rank, A is unique although V and W may 
not be. A is the parallel transporter we assign to each 
two-qubit link in our gauge theory. We give more details 
about the properties of S and A and their calculation in 
sections III CI and III Dl respectively. 

Just as one can also decompose an arbitrary matrix 
into the right polar form using the standard SVD, one can 
also do the same for the Lorentz SVD. That is, we may 
write S{a, b) = S{a, b)'A{a, b)' where S{a, b)' = VaSabVj 
and A(a, &)' = rjVariWff . The interpretation of this de- 
composition is that now A(a, b)' represents the local op- 
eration performed on qubit b that symmetrizes the cor- 
relations in the two qubit state. S{a, b)' is the (different) 
symmetrized state. While A(a, b)' ^ A(a, b) (they are re- 
lated by A(a, b)' = 'i]A{a, b)rj) the eigenvalues and trace of 
a cyclic product of the A representing the overall trans- 
formation around a loop remain the same provided one 
is consistent in taking just the left or the right polar de- 
composition. It is this total transformation that we look 
at in the next section. 

One might worry that restricting the matrices V and 
W to be elements of the Lorentz group, one may not 
always be able to bring a two qubit correlation ma- 
trix to Bell diagonal form and indeed some cases do 
exist. These non-diagonalizable cases are not generic 
however. Verstraete et al. have listed the non- 

diagonalizable cases, most of which turn out to be prod- 
uct states comprised of one or more pure subsystems. 
We do not consider two qubit links composed of prod- 
uct states since it seems unlikely one can consistently 
define a correlation related parallel transporter in such 
examples. There is one other non-diagonalizable case, 
it is the state po|00)(00| + pi\^+){-^+\ +p2|*")(*"|- 
These states can still be brought to a well defined diago- 
nal form with infinite Lorentz transformations, a process 
known as quasi-distillation, and we use this technique in 
section IIVI However, one does not have to apply infi- 
nite Lorentz transformations to assign a parallel trans- 
porter to these states; the parallel transporter only has 
to symmetrize {S = S'^) the correlations in the two qubit 
link. This can be done without the use of infinite Lorentz 
transformations. 



B. Twist as a measure of asymmetry of 
correlations: A Lorentz group gauge theory 

We now take a cyclic product of the parallel trans- 
porters A representing a loop around the lattice points 
(qubits) to form our measure of the degree of asymmetry 
of bipartite correlations globally. We define this measure, 
twist ^ as 

C(a6 • • • z) = 3tr{A(a, z) ■ ■ ■ A(c, b)A{b, a)}. (9) 

The constant 1/4 is chosen to give ^ = 1 when the over- 
all transformation is the identity. We now ask how the 
parallel transporters transform under local operations or 
in the language of gauge theories, gauge transformations 
U e SO+{l,3). From Eqs. @ and © we see that 
Va UaVa and WbUb- This implies A(a, b) 

transforms as 

K{a,b)^UaA{a,b)U-\ (10) 

just as parallel transporters should indeed transform un- 
der gauge transformations. We have used the identi- 
ties r]r] — 1 and rjU^rj = U^^. One can now see 
under gauge transformations, the parallel transporters 
do change but the total transformation around a loop 
only changes up to a similarity transformation leaving 
the trace and the eigenvalues of the overall transforma- 
tion A(a, z) • • • A(c, 6)A(&, a) invariant under 50+(l,3). 
This overall transformation is sometimes also known as 
a holononiy. In the case of twist it is a Lorentz group 
holonomy. 

Twist is therefore a measure of the asymmetry of the 
correlations you cannot 'gauge' away glob ally. One is 
reminded of the Aharonov-Bohm effect [21| and the geo- 
metric phase [31 as illustrations of features that although 
can be gauged away locally, cannot be gauged away glob- 
ally. Twist shares another property of the Aharonov- 
Bohm effect; that one cannot associate a twist to an in- 
dividual link since each parallel transporter associated to 
that link can be 'gauged away'. We can only meaning- 
fully associate a twist globally to the entire loop of links. 
In this sense, twist is a nonlocal property of the bipar- 
tite correlations present in the loop. We give a simple 
illustration of the idea in Fig. [T] 

Experimentally determining the amount of twist in 
an unknown quantum state is a similar problem to de- 
terming the amount of entanglement. Both can be inter- 
ferred from state tomography. If however, one knows the 
state beforehand, one can confirm the amount of twist 
by applying sequentially local operations to each qubit 
around a particular loop that symmetrize or 'untwist' the 
bipartite links. Using our three qubit example illustrated 
in figure [U one applies the local operation A(6, a)^^ 
to qubit b to symmetrize the link ba and then applies 
K{b,a)~^ K{c,b)~^ to qubit c to symmetrize be. One 
should find that the last link ac can only be symmetrized 
by the overall transformation A(6, a)~^A(c, 6)~^A(a, c)~^ 
on qubit a but in doing so makes the neighboring link. 
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FIG. 1: Graphical representation of twist ^(abc) for a state of 
three qubits. The ribbons between the qubits represent the 
'twist', the degree of asymmetry in the correlations in each 
two-qubit link. The twist in each link is not invariant by 
itself however. One can untwist each link by careful choice 
of local gauge transformation (representing a local operation) 
U £ SO"'' (1,3). Fig. |l(b) corresponds to the local transfor- 
mation on b that results in the link A(6, a) being the identity 
A{b,a) UbA{b,a). Similarly fig. 



1 c 



corresponds to the 

local transformation on qubit c that untwists A{c,b)Uj^^ i.e. 



A(c,fo)f/,- 



(7cA(c,b)C/-\ One can still perform a local 



transformation on qubit a however the total twist in the link 
will not change. Even though locally we can untwist each link 
we cannot gauge away the twist globally. 



a&, assymetric by the same amount. This could be con- 
firmed by tomography on one of the bipartite hnks. One 
should note since these operations are generally proba- 
bilistic they are not always successful. If one obtains the 
undesired outcome one would need to apply a compen- 
sating operation that cancels out this undesired outcome 
which is also generally probabilistic. 



C. Properties of A and S 

In the Lorentz SVD given in eq. (O we have chosen 
to make the three singular values si, S2 and S3 take the 
same sign as this results in a unique A when the singu- 
lar values are all non-zero. Giving them the same sign 
also seems natural as the three spatial correlation sin- 
gular values are treated equally. There are two possible 
cases: (i) The sign is negative if det5 < and (ii) the 
sign is positive if det S > 0. For case (i), the diago- 
nal matrix E represents a Bell state mixture, the largest 
component of this mixture being the singlet, l^*"). In 



case (ii), the singlet makes up the smallest component 
in the Bell diagonal mixture. In both cases the order- 
ing of the remaining three Bell states, |$^), |$^) and 
1 5*+) does not affect the form of the parallel transporters 
A however it does affect the form of the symmetrized 
correlation matrix S. One can see why the ordering of 
the singlet is the determining factor on the form of A 
from its special properties, namely that it is invariant 
under symmetric local operations i.e. A^A\'if'~) — \'i>^) 
where A G SL{2, C). This last fact can be seen even more 
clearly in the correlation matrix form since S'|,p-^ = 77 and 
our relation becomes equivalent to the defining property 
of the Lorentz group, eq. (|4|). That is, under symmetric 
local operations, the singlet state is invariant, it remains 
in the anti-symmetric state space, a state space consist- 
ing of a single point. To change the ordering of the singlet 
state in the Bell mixture one has to apply different local 
operations to the two qubits. The three remaining triplet 
Bell states are not invariant under symmetric, local op- 
erations, however they do of course remain in the much 
larger symmetric state space. In fact one can transform 
any triplet Bell state into any other using a subset of 
these symmetric local operations. Using these facts we 
can see that the ordering of the triplet states in a Bell 
state mixture does not affect the form of A since making 
the replacements V — >■ VU and W — WU leave it invari- 
ant. U e 50+(l,3) is a local operation applied to both 
subsystems. One also notes that the absolute singlet con- 
tent in S is unaffected by symmetric local operations. It 
is thus the ordering of the singlet that determines the 
form of the parallel transporters and the absolute singlet 
content in the symmetrized states. 



The parallel transporters are well defined and unique 
when S is full rank however they are not unique when one 
or more Si are zero. It is necessary but not sufhcient for 
S to be full rank for entanglement to exist in that link. 
Therefore the assignment of a unique parallel transporter 
to a link does not mean that link is entangled, however 
it does imply that the correlations in that link are 'quan- 
tum'. That is, correlations that do not appear in classi- 
cal physics. This type of weaker correlation is generally 
known in the literature under the name discord p2| . In 
this paper we only consider the cases where S is full rank 
and therefore the parallel transporters are unique. How- 
ever one may be able to define the notion of a partial 
holonomy in the cases of rank(5') > 2 corresponding to 
classically correlated (rank 2) and partially quantum cor- 
related (rank 3), by assigning a parallel transporter over 
the subspace in which the singular values are non-zero. 
This idea has been proposed in the context of geometric 
phases by Kult et al. [2j|. In the case S is rank 1, the 
bipartite link is a product state meaning that the two 
qubits are uncorrelated. The two qubits know nothing 
about each other and never can do if we restrict to using 
only local operations. Assigning a parallel transporter in 
this case does not seem possible. 
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D. Calculation of the parallel transporters, A 

Calculation of a parallel transporter for a given link 
involves finding the two elements, V and W in 5*0+ (1, 3) 
that diagonalize 5 in eq. ©■ A typical S will have all 
non-zero elements. An arbitrary element of 5*0"^ (1,3) 
may be decomposed into a product of two different el- 
ements of i?i,i?2 G SO{3), two spatial rotations, sand- 
wiching an element of i? e 5*0(1,1), a Lorentz boost 
along one spatial axis, V = RiBR2. Explicitly the 4x4 
representations are: 



B 



'I 

CqC/3 s 

CaS/3 S 

\Q -s„ 











coshiy9 sinhiy9N 

10 

1 

sinh(/3 coshc^y 



(11) 



(12) 



where = cos a and Sq, = sin a and analogous relations 
for the two other Euler angles (3 and 7. This decom- 
position suggests an iterative procedure to diagonalize S 
by applying a rotation R followed by a boost B to each 
qubit in the link. The aim is to depolarize each of the 
individual qubit 's Bloch vectors simultaneously making 
each individual qubit 's density matrix proportional to 
the identity. Rotations do not change the length of these 
Bloch vectors but they can change each of the spatial 
components. One strategy one might try is to first ro- 
tate each qubit 's Bloch vector so that it is aligned along 
one spatial axis, say the z axis, making each of the en- 
tries ^oi, S'02 and ^lo, 5*20 zero and then boost along 
the z axis in the opposite direction to which that Bloch 
vector points by a small amount. One repeats this until 
both qubits are simultaneously depolarized. The state 
is then proportional to a Bell diagonal state, although S 
may not be diagonal or in the correct signature. One can 
diagonalize the final 3x3 block by applying two final spa- 
tial rotations, More detailed calculational techniques 
are reported in [H [H Hi HI . 

Another strategy rather than finding V and W directly 
is to simplify the form of the two qubit links by apply- 
ing gauge transformations to each of the qubits since the 
overall transformation around a loop is gauge invariant. 
For example in the three qubit case, one may simplify 
calculation by choosing Ua, Ub and Uc to reduce the com- 
plexity of diagonalizing S i.e. make the transformations 



S{b,a) -> UbS{b,a)Ul, 
S{c,b) -> U,S{c,b)U^, 
S{a,c) -> UaS{a,c)U^. 



(13) 



One may also find it simpler to work in the density 
matrix basis; rather than dealing with real 4x4 matrices, 
one can find 2x2 complex matrices, A,Bg SL{2,C), 
that take each two qubit density matrix, p, to a Bell 
diagonal state, p'. That is p' = A<Si BpA^ (g) B'^ . 



Analytical calculation of V and W is possible when 
the singular values are non-degenerate corresponding to 
unique V and W. One can solve the eigenproblem for 
S^-qSr] = WY?W~^ . is a Lorentz matrix whose 

columns are the eigenvectors of W . Likewise one can 
find V by solving the eigenproblem for SrjS^rj. 



III. RESULTS FOR TWO-QUBIT STATES 

A. Pure states of two qubits 

We first consider a pure state of two qubits which can 
be written (up to local unitary equivalence) as a|00) -I- 
/3|11) where a and (3 are real. We wish to calculate twist 
for this state, the path from qubit a to 6 and back again 
given by 



^(a6) = itr{A(a,6)A(6,a)}. 



(14) 



We can use gauge transformations Ua and Uh to simplify 
the calculation of A. Using the SO^{l, 3) transformation 
on qubit a 



Ua 



( 



1 








2afi 





-1 














1 





2a0 








^2^1 



(15) 



and the transformation [/;, = I one can convert each link 
to S{a, b) = Eab = a/? X diag{l, —1, —1, —1}, a state pro- 
portional to the singlet, which, of course, has symmetric 
correlations under interchange of the qubits. This is the 
well known result that given a single copy of a pure, en- 
tangled state of two qubits, one can always convert it to 
a singlet using SLOCC [7]. The gauge transformation Ua 
is equivalent to the SL(2,C) matrix 







A = 







(16) 



in the density matrix basis. 

Since we can gauge away all the freedom in each link, 
the parallel transporters can all be consistently gauge 
transformed to A{a,b) = A(6, a)"-'^ = I giving £,{ab) = 1. 
This gives the result that all pure states of two qubits 
are untwisted. 



B. Mixed states of two qubits 

Arbitrary two qubit states are also untwisted (^ (ab) — 
1). This is simply seen by noting that S{a, b) = S{b, a)^ . 
Since S{a,b) = Val^abWif , S{b,a) = WbE.bVj. Note 
S(jb is now generally no longer proportional to the singlet 
but it does have the same signature (+,—,—,—). The 
decomposition of S{a,b) implies A(a, 6) ~ A{b,a)^^ or 
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equivalently there exist gauge transformations, Ua and 
Ut that untwist the hnk. For example, choose Ua — V^^^ 
and Ub — Wf^^. This illustrates another key property of 
gauge theories; one can only get the possibility of a non- 
trivial holonomy from a loop that encloses some area. 
Clearly for two-qubit paths this can never be done. One 
needs to consider states of three qubits and higher for 
the possibility to see non-trivial twist. 



IV. RESULTS FOR THREE QUBIT STATES 



A. Pure states of three qubits 



It is well known that a pure state of three qubits can 
be brought to one of two standard forms using SLOCC: 
(i) The GHZ state |000)-h|lll) or (ii) the W state |001) + 
|010) -I- 1 100). A generic pure three qubit state is of the 
GHZ class [5| . Since both of these states are symmetric 
under permutations of the qubits we expect that pure 
states of three qubits have trivial twist. However, we 
prove that the overall transformation is a tt rotation. 

Theorem: A generic state of three qubits in a 
pure state has $^(abc) = 0, with total transformation 
A(a, c)A(c, 6)A(6, a) = diai/ll, 1, — 1, — 1} in some basis. 
That is, the twist around the three qubit loop is a tt 
(Pauli) rotation around some spatial axis. 

Proof: We prove this by again making a careful choice 
of gauge transformation, C/, on each qubit locally. After 
this choice one finds that there is a Pauli rotation that 
cannot be gauged away globally. 

We use parameterizations of GHZ and W states first 
given by Diir et al. f| and solve for each class individu- 
ally. These two parameterizations populate the whole of 
pure three qubit state space. 

Any state in the SLOCC GHZ class can be written as 



\GHZ) = cs\mO) + sse'P\c^aMc 



(17) 



up to local unitaries. In the GHZ form 



|1), |0h) =c/3|0)-t-s^|l) and|0e 



,|o)- 



a) = Ca|0) + 

s.y|l). c and s 



stand for cos and sin respectively and the real angles have 
the ranges a, /3, 7 e (0, 7r/2], 6 G (0, it /A] and ip e (0, 27r]. 

We next calculate S{b,a), S{c,b) and S{a,c) for each 
two qubit link and try to get each link into the S form 
using local operations (gauge transformations) U. We 



choose the following SL{2, C) transformations 
/ 



A = 



B 



C = 



fc)^ 



piip „2 



(18) 



(19) 



(20) 



t here stands for tan and the corresponding Lorentz ele- 
ment Ua is given either by Eq. ^ or by 



Ua = T{A®A*)T'^ 
and likewise for Ub and Uc- T is the matrix 



1 



a 
1 1 



T = — 

a 



(21) 



(22) 



Using the gauge transformations Ua, Ub and Uc on 
S{b,a), S{c,b) and S{a,c) i.e. 

Sib,a)^UbSib,a)Uj, 
S{c,b)^UcS{c,b)U^, 
S{a,c) UaS{a,c)Ul 



(23) 



we can reduce each S to diagonal form. However, 
the gauge transformed UaS{a,c)U]: has the signature 
(+,+,—,+). We therefore need to multiply it by 
diag{l, — 1, 1, — 1}, a Pauli y rotation, to put the final 
link in the desired form. This final Pauli rotation is the 
holonomy, since we cannot gauge it away. In other words 



A(a, c)A(c, 6)A(6, a) = diag{l, -1, 1,-1} 



(24) 



giving ^{abc) — for all states in the GHZ SLOCC class. 

Note that the result would have been the same without 
the use of the gauge transformations since the overall 
transformation is gauge invariant. The careful choice of 
gauge transformation just made the calculation simpler. 

Any state in the SLOCC W class can be written as 



\W) = w|000) -I- x\001) + y|010) + z|100) 



(25) 



up to local unitaries with w, x,y, z > 0. We can follow the 
same strategy as for the GHZ case, trying to reduce the 
complexity of finding A by gauging away as much of the 
twist as possible. If we choose the gauge transformations 



A 

B 
C 



Q 



^ 





t. / - 
y J 



(26) 

(27) 
(28) 
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we can reduce S{b,a) and S{c,b) to the desired S form. 
S{a,c) is also diagonal but with the wrong signature 
(+,+,+,—). Again, we have given these transformations 
in SL{2, C) form. However, it is simple to find Ua, Ub and 
Uc in SO+{l,3) form using cither Eq. © or Eq. A 
key step in the gauge transformation procedure is using 
the local operator Q = diagjri, l/n}. To get these states 
in the desired Bell diagonal form we must take the limit 
n — cx) which takes Q to the limit of a projective mea- 
surement. Since projections are not invertible they are 
not group elements and strictly do not belong to repre- 
sentations of S'L(2,C) or S'0+(l,3). However, one can 
get infinitesimally close to a projective measurement and 
the Lorentz singular values of S are well defined in this 
limit. In fact each S becomes proportional to a pure 
Bell state, the links S{b,a) -^yzx diag{l, — 1, — 1, — 1} 
and S{c,b) — > x diag{l, — 1, — 1, — 1} being propor- 
tional to the singlet, l^*") and the final link S{a,c) — > 
xz X diag{l, 1, 1, — 1} being proportional to |^+). This 
is a process known as 'quasi-distillation' Two qubit 
states of this form can be brought to pure Bell states but 
the probability of success decreases to zero. 

These gauge transformations untwist all the links ex- 
cept for S{a, c) which has the incorrect signature. There- 
fore the overall holonomy for W states is again a Pauli 
rotation, A(a, c)A(c, 6)A(&, a) = diag{l, — 1, — 1, 1} {z ro- 
tation in this particular basis) giving S,{abc) = once 
again. 

Discussion: We find that the total transformation 
around a loop for pure states of three qubits is equiva- 
lent to the matrix diag{l, — 1, 1, — 1}. In the GHZ case 
we find that the links ab and be could be brought to 
the form Ao|*")(*"| + Ai|*+)(*+| where Aq > Ai by 
the local gauge transformations U whereas the remaining 
link ac can only be brought to the form Ao| <&"'") ('I'^l + 
Ai|<I>^)(<I>^|. The links ab and be are in the desired form 
i.e. they correspond to E with signature (-1-,— ,—,—), 
these are Bell diagonal states whose largest contribution 
in the Bell state mixture is the singlet. We have chosen 
this particular signature as it makes the parallel trans- 
porter A unique provided S is full rank and seems most 
natural (the spatial components of S corresponding to 
ax, <Jy and Oz are all treated equally). The singlet l^*^) 
is also unique amongst the Bell states since it is the only 
anti-symmetric state whereas the other three Bell states 
are symmetric. We found that one cannot put every two 
qubit link in the form Ao|*~)(*~| -|- Ai|*+)(*+| consis- 
tently with local gauge transformations. For the GHZ 
case we find two of the three links can be brought to 
mixtures composed of the singlet while the remaining 
link must necessarily be a mixture of the symmetric Bell 
states. One can see a similar thing happens with the W 
state case. 

Note that each of the two qubit density matrices have a 
maximal rank of two as a consequence of the Schmidt de- 
composition of the pure three qubit state. This fact gives 
one enhanced freedom in choosing local gauge transfor- 
mations U and allows one not only to symmetrize each 



two qubit state simultaneously but also put each link in 
Bell diagonal form maximizing the entanglement in each 
link. Rank two density matrices correspond to a S with 
form 

/y \ 
V _ I -X 

■ ^^^> 

Vo -y) 

The Lorentz singular values are given by y = 5(^0 + Ai) 
and X — |(Ao — Ai) where Ai are the eigenvalues (given in 

non-increasing order) of ^ ® cfy)p* (cy ® (Jy)^, a 
familiar operator used in calculation of the concurrence 
Q. Now one can see where the extra gauge freedom 
comes in. E is composed of two degenerate blocks, de- 
generate with regard to the Minkowski metric. In other 
words S is the direct product ycr^ © ~xl. This means 
E is invariant under transformations of the form uYjV?^ 
where u = 5*0(1, 1) ® 5*0(2) giving the extra freedom in 
choosing the last gauge transformation. 

The interpretation of the result prompts the question 
whether one can put each of the three two qubit links 
in a mixture of symmetric Bell states simultaneously. A 
mixture of symmetric Bell states for the case of rank two 
Bell diagonal states corresponds to E = diag{y, y, x, — x} 
{p = Ao|^'+)(«'+| + Ai !$+)($+ 1). One can indeed find U 
that take each link to a state of this form simultaneously. 
However, it should be noted that the parallel transporter 
assigned to each link A — VrjW"'"r] is no longer unique 
since E no longer has the Minkowski signature. We can 
however choose the A closest to the identity, in some sense 
the 'most parallel' transformation between the qubits. 
Including this extra rule one finds the overall transfor- 
mation associated to this new rule of parallel transport 
is the identity. 

As a final comment, we note that our result has similar- 
ity to the Wigner angle or Wigner rotation [23, . The 
product of two pure Lorentz boosts in different spatial 
directions is not equal to another pure Lorentz boost, 
rather it is a pure boost multiplied by a rotation, the 
Wigner rotation. Aravind has shown that the Wigner ro- 
tation can be thought of as a holonomy in rapidity space 
(hyperbolic space) resulting from the area enclosed by 
the two pure boosts and the pure boost connecting the 
end points. The magnitude of this rotation is given by 
the area of a hyperbolic triangle, the length of whose 
sides are given by the rapidities of each of the boosts. 



B. Mixed states of three qubits 

Mixed states of three qubits show a wide range of pos- 
sible holonomies. Generically, they have a full 50''"( 1,3) 
holonomy group structure however they may be hard to 
solve analytically. Here we give some simple examples. 
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1. Untwisted states 
There is an untwisted state 

\ (I* ® Ic + I* ® la + |*ac>(*acl ® ^b) ■ 

(30) 

For this particular state each two qubit hnk is the Werner 
state 

il+i|vI/->(*-|= (31) 

\\^-){^-\ + i (|*+>(vi/+| + |$+)($+| + |a>-)(<f>-|) , 

corresponding to the correlation matrix 5* = ^ x 
diag{l, — |, — ^, — ^}. Each of the links has no entan- 
glement, although this state is right on the cusp of being 
entangled - adding a infinitesimal amount more of the 
anti-symmetric Bell state l^*^) to the mixture would re- 
sult in an entangled state. One also notes this state is 
an equal mixture of symmetric and anti-symmetric Bell 
states. 

There is also the complementary untwisted state 

^{}GHZ){GHZ\ + \W){W\) (32) 

where the states are \GHZ) ^ y^|000) -I- ^§1111) and 

\W) = y^(|001) + |010) + |100)). Each two qubit cor- 
relation matrix has the form S — ^ x diag{ 1 , ^ , ^ , ^ } 
i.e. it is the same as the example above except each 
link is in the signature (-I-, -|-, +), that is each link 
is an equal mixture of the three symmetric Bell states, 
|$+)(<I>+| + |<I>-)(<I>-| + |*+)(*+|. This state is conjec- 
tured to be the state of two qubits with maximal disso- 
nance, that is, maximal quantum correlations that are 
not entanglement (29j . 



2. 5*0(1, 1) holonomy from rank 3 bipartite density 
matrices 



The three qubit state 



p\GHZ){GHZ\ + {l-p)\^){i,\ (33) 



shows a wide range of interesting behavior depending on 
the values of its variables. Depending on these values, 
the overall transformation around the loop can either be 
the identity, a tt rotation or a general element of 5*0(1, 1). 
Each link can exist in three possible regions in parame- 
ter space, two having entanglement in the two qubit link 
and one region being separable. The regions are sepa- 
rated by two critical points such that one of the singu- 
lar values Si = resulting in det S ^ 0. The critical 
points in this example mark the transition from the link 
being entangled to separable. The pure states are de- 
fined in eq. ^ are \GHZ) ^ ^TjdOOO) -I- |111)) and 

IV") = 2;|001) + y|010) -I- z|100) + w|lll). We choose the 
amplitudes of \ip) to be real and positive. 

We solve this example analytically by finding the form 
of the diagonalizing matrices V and W for each link di- 
rectly without prior use of gauge transformations. One 
may then have to apply a Pauli rotation to get E into the 
correct signature. The Pauli rotation applied depends on 
which of the regions the bipartite link belongs which in 
turn depends on values of the variables p, x, y, z and w. 
We solve for the link a6, solutions for the other two links 
follow in an analogous way. The form of the correlation 
matrix is 



S{h, a) 








\(1 -p)i-w'^+x'^ 





2(1 — p){wx + yz) 





{1 ~ p){-w'^ + x'^ + y"^ ~ z"^) 


—2{l—p){wx^yz) 

p+{l-p){w'^ + x'^ -y'^ - z"^)/ 

(34) 



One may diagonalize this correlation matrix with Vb ~ 
T{B ® B*)T'^ and Wa = T[A ® A*)T^ where 

[ mn 

and e{q) = \/ p + 2(1 — p)q^. The resulting diagonal ma- 



trix, S, has singular values 

So = ^e{w)e{x) + {l-p)yz, 

si = {1 ~ p){yz + wx), 

52 = (1 -'P){yz - wx), 

53 = ^£{w)e{x) ~ {I ~ p)yz. (36) 

However, depending on the values of the variables the 
signs of the singular values may not be correct. There 
are three regions, /, // and ///, corresponding to the 
different signs S2 and S3 may take, sq > si > for all 
values. Region / corresponds to the range of parameters 
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in which wx > yz, with the first critical point occurring 
at wx = yz when S2 — 0. This region is entangled with 
conciuTcnce C/ = 2(1 — p)(wx — yz) = —2si however 
the signature of S is (+,+,—,+). We can bring it to 
the correct signature (+,—,—,—) by applying a Pauli y 
rotation to either A or B, this choice does not change the 
form of the parallel transporter A we assign to the link. 
If we choose to put the rotation on qubit o we replace 
A iAoy . 

Region // is separable although it still has dissonance 
and this occurs over the range wx < yz < ^2(i-p^^ cor- 
responding to dctS > 0. In this case the signature of S 
is correct, (+,+,+,+), and we can construct the paral- 
lel transporter from V and W without applying a Pauli 
rotation. 

Region III is also entangled and lies in the range yz > 
'^2('i-p) ■ '^^^ concurrence in this region is given by Cm = 
2{l—p)yz—e{w)e{x) = — 2S3. The signature of this region 
is {+, +, +, — ) therefore one needs to apply a Pauli z 
rotation to one of the qubits. 

The parallel transporters one assigns to each link are 
thus determined by the region that each two qubit link 
resides in. There are 3^ = 27 possible combinations as- 
signing three regions to three links and therefore 27 pos- 
sible overall transformations (or holonomies) around the 
three qubit loop. Many of these combinations turn out 
give the same overall transformation however. 

We investigate the possible forms by first making some 
observations about V and W. They consist of two blocks, 
one spanning the identity and z component which be- 
longs to the group 5*0(1, 1) and a second block equal to 
the identity spanning x and y. If the coefficients in the 



link are such that it belongs to region / we need to ap- 
ply a Pauli y rotation to get the signature correct. This 
rotation takes one oiV or W to the block diagonal form 
0(1, 1) ® CT^ and thus the parallel transporter assigned to 
a link residing in region I also has this form. Applying 
the same arguments to regions II and ///, we find 

Ai = 0(l,l)©a„ 
An = SO(l,l)eI, 
Am = SO{l,l)®-I. (37) 

We can therefore observe that any loop with an odd num- 
ber of links in the entangled region / will have an overall 
transformation belonging to 0(1, 1) © cr^. An element 
from this group can only have eigenvalues {1,-1,1,-1} 
giving ^{abc) = 0. That is, the overall transformation is 
a TT rotation. 

Calculation for the case when all three links are separa- 
ble (all three in region II) yields a trivial total transfor- 
mation, the identity, corresponding to £,{abc) = 1. Prom 
this case one can infer that combinations of region II and 
region III result in a total transformation of either the 
identity {^{abc) = 1) or a tt rotation {^{abc) = 0), the 
former case occuring if there is an even number of region 
/// links and the latter if there are an odd number of 
region /// links. 

The only remaining cases are those that feature an even 
number of region / links. These cases have non-trivial 
holonomy. Wc first look at the cases where all links are 
entangled that is A/(a, c)Aj(c, 6)A///(6, o) and the two 
other permutations of region III. The eigenvalues of 
this transformation are 



A/ (a, c)Aiic, b)Am{b, a) = diag {ff{f{, -1,-1, 
AKa, c)Ajjj{c, b)Ai{b, a) = diag {gg, -1, -1, 
Am{a, c)A/(c, 6)A/(6, a) = diag {f^, -1, -1, 



} ^ = ^^f^ffe^, (38) 

} e(a6c)= <-^^i^;;gy , (39) 

} ^i-bc) = ^ttaZT ■ (40) 



The final case, that of one unentangled link, 
A/(a, c)Aj{c. b)Aii{b, a), and its permutations arc identi- 
cal to the case above, however the eigenvalues taking the 
value —1 become 1. 

In this example the holonomy is determined by the lo- 
cation of the crictical points in each link. These critical 
points mark the transition from entangled to separable. 
One can verify analytically that any two qubit density 
matrix with rank three or smaller will also have critical 
points determined by the entangled to separable tran- 
sition. However, for two qubit density matrices of rank 
four this behaviour no longer holds as we illustrate in the 
next example. 



The special role of the singlet in the holonomy is also 
highlighted again in this example. It is interesting and 
strange to note one can put all links simultaneously into 
states that do not have any singlet content. That is, one 
can make each S a mixture of the three symmetric Bell 
states Ao|*+)(\E'+K Ai|$-)($-| A2|$+)(*+| and the 
holonomy will be the identity no matter which of the 
regions each of the two qubit links resides. This is not 
possible if one requires any singlet content, the holon- 
omy is always a tt rotation or a more general element of 
50(1,1). We believe this is due to the malleability, or 
size of the the symmetric state space; the ordering of the 
the three triplet states does not change the form of the 
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parallel transporters. The singlet is much more rigid, its 
ordering will change the parallel transporters. 



SO{l, 1) holonomy from rank 4 bipartite density 
matrices 



An example of a state with S'0( 1,1) holonomy with 
two qubit density matrices that are each full rank is the 
following state 



p\W){W\ + [1 - p)\W){W\. 



(41) 



The state \W) is defined in eq. ([25]) and \W) = crf^|W^). 
In the following we take w = (choosing w results 
in a SO+{l, 2) holonomy) and p g [0, 1]. 

For these states we chose to solve for the parallel trans- 
porters directly without the prior use of gauge transfor- 
mations. That is, we find the V and W for each two qubit 
link in the Lorentz SVD given by eq. ([6]). One finds the 
subspace of each S spanned by ax and ay is already diag- 
onal and the singular values associated to this subspace 
are degenerate. Only the subspace spanned by I and a^ 
needs to be diagonalized. This can be achieved by simple 
Lorentz boosts applied in the z axis given by the group 
SO{l, 1). Explicitly each V and W has the form 



/coshiy9 sinhiysN 

±1 

±1 

\sinh(^ coshipj 



(42) 



and therefore the parallel transporters and the overall 
transformation around the loop have the same form. 
That is, the holonomy is a simple pure boost along some 
spatial axis. Two of the singular values for each link 
are degenerate and take the same sign. The last spa- 
tial singular value associated to the boost however can 
change sign. It is negative when det 5 < in which 
case the central diagonal block in eq. (|42|) becomes —I 
and I when det S > 0. This critical point occurs at 
x^y^p{l — p) = \/p{l — p){y'^ — + y^z"^ for the link 
ab with analoguous relations for the other two links. 
There are therefore two regions seperated by the sign 
change of this singular value. The sign change in this 
example does not mark the transition between entan- 
gled and separable since the two qubit density matri- 
ces are rank 4 rather than rank 3. However, we can 
relate this point to (negative) concurrence; at the critical 
value the Bell state with the smallest probability in mix- 
ture represented by E has probability p^ = — 2C where 
C = Ao — Ai — A2 — A3 is negative at the critical point. 

It is not hard to solve for this case. One finds the twist 
is given by hyperbolic cosine of the sum of the hyperbolic 
angles, ip, 



^{abc) — cosh^ 



2 / fac + Vcb + fba 



when there are an even number of links with det S < 
and 



^{abc) = sinh^ 



2 / fac + fcb + (Pba 



(44) 



when there are an odd number of links with det S* < 0. 
The rapidities (hyperbolic angles), ip, familiar variables 
from special relativity, are given by 



1. fl + (3 



and each /3 is 



Jba 



Pcb = 
Pac = 



(l-2p)(y2-z2) 

{l~2p){x^~y^) 
{l-2p){z'' -X^) 



(45) 

(46) 
(47) 
(48) 



(43) 



Notice when p g {0, i, 1}, ^ = 1 i.e. one can untwist the 
bipartite correlations globally. One can also untwist the 
bipartite correlations for all p if two or more oi x, y and 
z are equal since the state before gauge transformations 
is more symmetrical. 



V. CONCLUSIONS 

In this paper we have introduced a new measure quan- 
tifying the global degree of asymmetry of two qubit cor- 
relations in the spirit of a lattice gauge field theory. We 
call this measure twist and it is phrased as the famil- 
iar Wilson loop. Twist is defined as the trace of the 
total transformation or holonomy around a loop, a fea- 
ture that cannot be 'gauged' away by general local op- 
erations, the gauge group naturally emerging being the 
Lorentz group 5*0^(1,3) familiar from special relativity. 
This measure can be interpreted operationally in terms 
of an SLOCC protocol. We defined this measure as the 
trace of the total transformation which is a multipartite 
5L(2,C) invariant of the state. However one may be 
able to get more information about the state by consid- 
ering the eigenvalues of the total transformation which 
are also multipartite SL{2, C) invariants. We point out 
that these invariants are not the polynomial invariants 
that have been typically studied in the literature. We 
have found that pure states of three qubits have a total 
transformation that is a tt rotation and provided some an- 
alytical examples of mixed three qubit states with more 
general transformations, elements of 5*0(1, 1), the group 
of Lorentz transformations in one dimension. 

In this work we have been investigating the possible 
analogy between nonlocality in well known gauge field 
theory effects such as the Aharonov-Bohm effect and the 
nonlocality present in certain quantum states violating 
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Bell inequalities. The aim being to phrase the nonlo- 
cality found in entangled states in terms of the familiar 
geometrical and topological ideas from gauge field the- 
ories. Our twist invariants are nonlocal, but not in the 
sense of states violating Bell inequalities. They are non- 
local in the sense that one may not attribute a twist to 
an individual link in the loop, since one may gauge it 
away to the trivial transformation. The twist is property 
of the entire loop of links and thus measures the global 
asymmetry of the two qubit correlations in that loop. 

We have found that entanglement is not necessary to 
be able to define a unique parallel transporter to a link 
except in the case of rank 2 and lower two qubit density 
matrices. In the rank 3 case, entanglement does however 
define the crictical points in the link and these points af- 
fect the parallel transporter one assigns to a link. Gener- 
ically however, to assign a unique parallel transporter to 
a link, one needs a weaker type of quantum correlation 
known in the literature as discord. Two qubit entan- 
glement is determined by the magnitudes of the dilated 
correlations. One can phrase it as the sum of the sin- 
gular values of the correlation matrix. Twist however, 
appears to characterize the asymmetry of the quantum 
correlations. In constructing twist one makes use of the 
additional information about the state encoded in the lo- 
cal operations that bring the correlation matrix to Bell 
diagonal form. 

The constructions used in this paper are all based on 
two qubit correlations. It is not clear that these ideas 
generalize simply to d level systems or three and higher 
qubit correlations. For a two qubit correlation matrix 
there is a simple standard form; one can bring any two 
qubit state to a Bell diagonal state with local operations. 
It is not clear that this standard form exists for two qudits 
or three or more qubits. In the case of three qubits one 
has a 4 X 4 X 4 correlation tensor. One can imagine several 
possible standard forms for such an object. 

We have several questions motivated by the similarity 
of this work to other ideas in modern theoretical physics. 



For example, in loop quantum gravity, the primitive vari- 
ables are also phrased as Wilson loops, however parallel 
transporters are assigned from the relation between the 
discrete points of space rather than the relation between 
qubits [30] . One possible avenue of investigation could be 
to look for states with bipartite correlations that mimic 
simple curved spacetimes. For example can one find a N 
qubit state that has the property such that the larger the 
loop around a central qubit becomes the total transfor- 
mation becomes more trivial? This could be a toy exam- 
ple of a simple curved spacetime such as the spherically 
symmetric Schwartzchild metric. 

One may also look for links between knots and quan- 
tum correlations. Witten has phrased Wilson loops in 
gauge field theories with a Chern- Simons action in 2 + 1 
dimensions as knot invariants [3l| . One wonders whether 
one can find states where a similar analogy can be made 
in our construction. 

On a more concrete note, we hope that the ideas pre- 
sented in this paper are useful in constructing and gain- 
ing an intuition for invariant correlation properties in N 
qubit states. 
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